Elastic domains in antiferromagnets 
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We consider periodic domain structures which appear due to the magnetoelastic interaction if the 
antiferromagnetic crystal is attached to an elastic substrate. The peculiar behavior of such structures 
in an external magnetic field is discussed. In particular, we find the magnetic field dependence of 
the equilibrium period and the concentrations of different domains. 



Spontaneous breaking of any discrete rotational sym- 
metry in monocrystals leads to the appearance of degen- 
erate phases characterized by different orientations of the 
order parameter. The phase boundaries have a positive 
(free) energy, otherwise the system would be thermody- 
namically unstable. Thus, in equilibrium one should ob- 
serve a single phase. This is, however, not the case for 
the phases with spontaneous magnetization (or electrical 
polarization). Periodic domain structures are formed in 
such materials in order to reduce the large, proportional 
to the volume, energy of the magnetic (or electric) field. 

Using the electron spin resonance technique, Janossy 
et al. 0,01 demonstrated the existence of the equilibrium 
domain structure in antiferromagnetic Y Ba2C'u30e+x 
and in Yi^xCaxBa2Cu30e. The domain structure was 
easily modified by the external magnetic field, and re- 
stored after switching the field off. The resonance method 
does not provide information about characteristic sizes 
and the arrangement of the domains. Recently Vinnikov, 
et al. j^l directly observed periodic domain structures 
in antiferromagnetic TbNi2B2C by the "finest magnetic- 
particle decoration technique" Q. Neither spontaneous 
magnetization, nor electrical polarization is seen in this 
material. Nevertheless, regular and reversible (by the 
change of the temperature and magnetic field) domain 
structures are realized with a periodicity of the order of 
a few microns in the samples of plate shape with a thick- 
ness 0.5 mm. 

In this Letter we discuss a possible scenario for the 
appearance of such periodic domain structures and the 
behavior of the structures in an external magnetic field. 
If a monocrystal is attached to some elastic substrate, the 
domain structure should inevitably arise in any orienta- 
tional phase transition in order to minimize the strain 
energy. Indeed, the stress free (but not strain free) mon- 
odomain state is realized if we discuss free surface bound- 
ary conditions. If the monodomain crystal is attached 
to some elastic substrate, stresses arise, and we loose a 
large, proportional to the crystal volume, energy. In this 
case, the appearance of the domain structure, with a pe- 
riod much smaller than the crystal size, will drastically 
diminish the energy of elastic deformations. 

All materials studied in the experiments of Ref. 0,0,01 
are collinear easy plane antiferromagnets with tetragonal 



FIG. 1: Schematic presentation of the domain structure in the 
antiferromagnetic crystal attached to an elastic substrate. 



crystal symmetry. For definiteness' sake we discuss this 
case. For the description of the antiferromagnetic states 
we use the exchange approximation 0|, treating the ef- 
fects of the magnetic anisotropy, magnetoelastic effects 
and the external magnetic field as perturbations. By a 
proper choice of the orientation of the coordinate system 
(say, the x axis can be oriented along the [110] or [100] 
direction) we can write the anisotropy energy in the easy 
xj/-plane as 



(1) 



and ly are the com- 



with a positive constant /3; here 
ponents of the antiferromagnetic unit vector 1. The 
four states l'^ = ly = 1/2, differ by the signs of the x,y- 
components (/*, \, \), and correspond to the min- 
imum of the energy, Eq. (Q. 

The elastic strain arises in the crystal due to the mag- 
netoelastic interaction 



lllxlyUxy + l2{l'l ~ ^Vii^xx 
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where Uik are the components of the strain tensor. For 
simplicity, we consider the magnetoelastic effects to be 
small compared to the anisotropy, Eq. The second 
term in Eq. ^ is zero for the discussed states with 
l'^ = ly. The elastic energy of tetragonal crystals can be 
written as a sum of 6 invariants: 

Pel Ml"ly + M2(uL + "yj + fJ'SiUxx - %y)^ + 



/i4ML + ^J,5'U'zz{Uxx + Uyy) + fle{Uxx + Uyyf- . (3) 

Minimizing the sum of the energies, Eq. (O and Eq. JSJ, 
we find that, in the monodomain state of the unstressed 
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FIG. 2: The displacement field inside the crystal far from the 
crystal-substrate boundary. 



crystal, the only nonzero component of the strain tensor 
is 
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Consider the antifcrromagnctic crystal, of plate shape, 
attached to a flat clastic substrate. Let the main axis 
of the crystal be normal to the plate. The domain 
structure that appears in this case is schematically pre- 
sented in Fig. 1. The state of each domain far from 
the crystal-substrate boundary is one of the unstressed 
ground states. Additional energetic contributions arise 
from i) the energy of the domain walls, and ii) the elastic 
energy localized near the crystal-substrate boundary in 
a layer of thickness of the order of the domain structure 
periodicity d . Minimization of these two contributions 
with respect to d leads to the equilibrium period of the 
domain structure. This argumentation is close in spirit 
to arguments of Ref. (see also |ill[l3). 

In accordance with the symmetry of the domain struc- 
ture. Fig. 1, the displacement vector u has only one com- 
ponent, Uy, inside the crystal, as well as in the substrate. 
Therefore, elasticity equation in the crystal reduces to 



Hid^Uy + ^i2c)luy = 0. 



(5) 



We assume that the substrate is elastically isotropic: 
^J'l = fJ'2 = 2//, where fj, is the shear modulus of the sub- 
strate. The boundary conditions at the sample-substrate 
contact surface are 

Uy\z =+0 = Wy|2 = -0,M29zUy|z = + = '^HdzUylz^^Q. (6) 

The displacement field inside the crystal, far from the 
crystal-substrate boundary, is a zigzag function u'-' as 
schematically presented in Fig. 2. It's expansion in 
Fourier series is 



0/ \ • 27r(2r7, + l)a; 

Tl = 



(7) 



with coefficients 



_ (-l)"7ia 



7r2^i(2n + 1) 



It is convenient to introduce a new field u"*" by the relation 
Uy — u° {x) + u^l^ {x , z) (8) 
The Fourier scries for in accordance with Eq. (j^J 



27r(2n+ l)x 



u+ ^ u+ sin 
The Fourier scries for u. 



■ exp 



27r(2n + l)y/jHz 



V 



in the substrate is 



27r(2n+ l)x ( 2TT(2n+l)z 
sm : exp 



n=0 



d { d 

From the boundary conditions, Eq. JHJ, we find 



2/i + y/jHjl2 



2m + \/7*i7^ 



(9) 



Now we can determine the total energy of the system 
per unit area in the xy-plane 



where ctq is the domain wall energy, and 



7r3 ^ (2n+ 1)3 

n=0 ^ ' 



0.27. 



The first term in Eq. (|10() represents the energy gain 
due to the arising of the domain structure. In the mon- 
odomain state of the sample with the free surface this 
is the only term of the magnetoelastic relaxation. The 
second term is the elastic energy due to the inhomoge- 
neous strain near the crystal-substrate boundary. The 
third term accounts for the domain walls energy. Mini- 
mization of the energy, Eq. (^UJ, with respect to d gives 
the equilibrium period 



d = dn — 
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(11) 



The classical law d oc i^/^, in the limit d <^ L, is the 
crucial point for the experimental verification of the do- 
main structure theory. In the case of antiferromagnetic 
domains there exists an additional valuable experimental 
possibility, namely to study the behavior of the structure 
in the external magnetic field. 

The magnetic field changes the orientation of the an- 
tiferromagnetic vector. The corresponding term in the 
energy is 



IHl 



(12) 



where X|| (Xi) is the magnetic susceptibility parallel 
(perpendicular) to the vector 1. In coUincar antiferro- 
magnets X|| < Xj.- 



FIG. 3: Two types of ordered domain structures, a) and b) 
and a disordered domain structure, c). 



FIG. 4: Dependence of tlie period of tfie domain structure on 
tfie magnetic field. 



If the magnetic field H is oriented along the y- 
direction, the minimization of the sum of the energies, 
Eq. d and Eq. gives ll - = /i^ for H < and 

ly^O for H > H,; here h = H/H,, = 2/?/(x± - X||). 
Firstly, the magnetic field leads to the appearance of ho- 
mogeneous deviatoric stresses a^x = ~^yy — 72 (^^ ~ 'y) 
due to the second term of the magnetoelastic energy, Eq. 
(O. These stresses are equal in the coexisting domains, 
and do not affect the periodicity. The second effect is the 
renormalization of the equilibrium strain tensor, Eq. 10}, 
by the factor (1 - /i4)i/2^ xhe third effect is the renor- 
malization of the energy of the domain walls. 

In order to find the domain wall structure one should 
take into account the energy of the spin space inhomoge- 
neous rotation 



(13) 



where 5 is a constant of the exchange interaction and 
is the rotation angle of the unit antiferromagnetic vector 
(Ix = cos(^). The variation of the sum of energies, Eqs. 
(|1I12I13() . gives the equilibrium equation 

= {h^ + siv? kp — cos^ ip) sin ip cos (14) 

where 5 = (g/2/3)^/^ is the effective thickness of the do- 
main wall. The first integral of this equation is 



sin^ tp 



(15) 



The constant of integration is defined by the condition 
far from the domain wall where dxf = 0. 
Finally, we find the domain wall energy 



'^±{h) = ctq I VT — /i^± 2/1^ arcsin 



(16) 



where (Tq = /3<5 is the energy of the domain wall without 
magnetic field. The minus sign in Eq. H16|) corresponds 



to the domain walls presented in Fig. 3a. These walls dis- 
appear at the critical field: ^ | ^ | ^ | ^ — s- — > ^ 
because the difference between neighboring domains van- 
ishes and the rotation in the wall tends to zero. The plus 
sign corresponds to the domain walls presented in Fig. 
3b. These walls remain at h > I: ^ | ^ | ^ | ^ . If the 
field is oriented along the x direction, the behavior of the 
domain walls in Fig. 3a, b will be opposite. The energy 
of the walls of type a) will be and at the critical field 
we obtain the following domain structures T I i I T I i • 
The walls of type b) will have energy cr~ , and vanish at 
the critical field t | T | T | T=^TTTT. 

In the general case one starts from an arbitrary distri- 
bution of the types of walls, as it is presented in Fig. 3c. 
It leads to a complicated and even irreversible behavior 
of the periodicity for increasing and decreasing magnetic 
fields. Nonetheless, a simple behavior can be achieved if 
one starts by applying a field higher than He, say in the 
y-direction. In this case, in equilibrium, one obtains a 
homogeneous state (say Ix = 1). By decreasing the field 
to H < He, a domain structure will appear with energet- 
ically preferable walls, with the minus sign in Eq. H16() 
[lll |. Then, the field dependence of the period is given by 



d-{h) 



h-^ - 2h^ arcsin 



1/2 



see the lower curve in Fig. 4. This behavior will be 
reversible in the magnetic field. If then, after turning off 
the field, one applies a field in x-direction, the domain 
walls will become unfavorable, and with increasing field 
the period will increase in accordance with 



d+{h) 



VT~~/i4 + 2h^ arcsin 



1/2 



see the upper curve in Fig. 4. But, if one stops the 
increase of the field at some value H < He, and starts to 
decrease it, the period will decrease more rapidly than the 



4 



upper curve, because it is favorable to produce new walls 
with the minimal energy. Note, that if magnetoelastic 
effects are comparable with the anisotropy, Eq. (Q, all 
obtained results for the behavior in the magnetic field 
remain valid. One merely should renormalize the value 
of He, and 6. 

The behavior of the domain structure will be much 
more complicated if one applies a field with some arbi- 
trary orientation. Then the nearest-neighbor domains 
should have different widths. Let us calculate the con- 
centrations of different domains, and c~ {c^ +c~ ~ 1), 
in the limit of small magnetic fields, Hx, Hy « He and, 
as before, assuming magnetoelastic effects to be small 
compared to the anisotropy. In this approximation we 
can neglect the small rotation of the domain structure, 
which inevitably arises in the general case, and also set 
= « 1/2. Since the average strain is zero, we have 
c~ = (i = x, ?/). At the domain 



Uyy = and u^iC^ + u 
wall we have boundary conditions of mechanical equilib- 
(T~j = axi and phase equilibrium |l2l | 



rmm, ct+ 



- axfrUi 



- 2(TxyUxy 



F- 



- Gxx^T 



'^'^xyUxy, (17) 



where F = F^e + F^i + Fh; <Jxx = dF/duxx and 
^xy = \dF / duxy are the components of the stress tensor. 
Solving this system of equations, we find the concentra- 
tions of the different domains: 



1 ± 



HxHy 

me 



,H, 



2 

me 



« m. 
(18) 

The domain period cannot be found analytically at arbi- 
trary field orientations. In this case there are no symme- 
try arguments forcing the domain walls, near the sample- 
substrate boundary (at distances of the order of the pe- 
riod), to be flat and oriented along the z-axis. The do- 
main walls should be inclined and curved in this case. 

Let us estimate the parameters of the theory. The 
anisotropy energy parameter is P^a^Ua"^, where 
a ~ 10~^ is the fine structure constant, U ~ lO^if is 
the atomic energy, and a ^ lO^^cm is the atomic size. 
The constant g ^ J /a, where J ^ IQ^K is the exchange 
energy. Then, for the domain wall width we find 5 ^ 



1/2 



aa 



cm. The shear modulus 
of a typical material is pi ^ Ua~^, and the magnetoelas- 
tic coupling constant 7 ~ a^Ua~^. Finally, for a crystal 
with thickness L ^ 0.5mm we obtain a reasonable esti- 
mation of the period d ^ \/6L to be of the order of a few 
microns. 

Note, that for easy plane antiferromagnets with hexag- 
onal crystal symmetry, anisotropy effects in the plane are 
small (^ a^) compared to magnetoelastic effects Eq. (0) 
(71 = 472). If we neglect this anisotropy the orientation 



of the domain structure will be arbitrary in the absence 
of a magnetic field. The antiferromagnetic vector will 
be oriented at an angle ±7r/4 to the domain walls. In 
the presence of the magnetic field the domain structure 
will presumably be oriented perpendicular to the field 
[ll| . The critical field is He ~ 7/ \/m(xj^ ^ X||) in this 
case. The domain concentration remains 1/2, and the 
dependence of the period on the magnetic field will be 
described by the value of d{h) obtained in the tetragonal 
case. 

This work is supported in part by the Deutsche 
Forschungsgemeinschaft (Grant No. DPP 1120), by the 
RFBR Grants No. 06-02-16509, 06-02-17281, and by the 
RF President Program. V.I.M. thanks Forschungszen- 
trum Jiilich for its hospitality. 



[1] 
[2] 

[3] 
[4] 

[5] 

[6] 



[7] 

[8] 

[9] 
[10] 

[11] 



[12] 



1176 (1999). 
Phys. Rev. Lett. 



91, 



A. Janossy et al., Phys. Rev. B 59 
A. Janossy, T. Feher, and A. Erb, 
177001 (2003). 

L.Ya. Vinnikov et al, JETP Lett. 77(9), 502 (2003). 
L.Ya. Vinnikov, LV. Grigor'eva, and L.A. Gurevich, 
Springer Ser. Mater. Sci. 23, 89 (1993). 
A.F. Andreev, and V.L Marchenko, Sov. Phys. Uspekhi 
23, 21 (1980). 

The additional elastic energy term appears only near the 
crystal-substrate boundary and does not appear near the 
free surface. This is a crucial difference between these 
elastic domains and the domains in ferromagnets. 

A. G. Khachaturyan, Theory of Structural Transforma- 
tion in Solids (Wiley, New York, 1983). 

B. Horovitz, G. R. Barsch, and J. A. Krumhansl, Phys. 
Rev B 43, 1021 (1991). 

B. Horovitz, Phys. Rev B 69, 140501 (R) (2004). 

E. Gomonay, and V. M. Loktev, J. Physics : Condens. 

Matter 14, 3959 (2002); Phys. Solid State 47, 1755 

(2005). 

The global domain wall orientations [x, or y) remain de- 
generate in the considered approximation. However, if we 
take into account demagnetization effects, the orientation 
of the domain structure with walls normal to the applied 
field will be more favorable. 

In the general case, the differential of a potential that 
is continuous across the interface contains only differen- 
tials of quantities which are also continuous across the 
interface: du-rr, dann, darr, dHr,dBn (n and r are the 
normal and tangential directions with respect to the in- 
terface and B is the magnetic induction). This statement 
is a generalization of the Privorotskii-Azbel condition for 
the coexistence of magnetic phases (see §44, in L.D. Lan- 
dau, E.M. Lifshitz, Electrodynamics of continuous media, 
Elsevier, NY, 1985) to the case where elastic effects are 
important. In Eq. 1171 we have neglected demagnetiza- 
tion effects (x << 1). 



